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Neutron scattering experiments have been performed on the ternary rare-earth diborocarbide 
Ce^^B2C2. The powder diffraction experiment confirms formation of a long-range magnetic order 
at Tn = 7.3 K, where a sinusoidally modulated structure is realized with the modulation vector 
q = [0.167(3), 0.167(3), 0.114(3)]. Inelastic excitation spectra in the paramagnetic phase comprise 
significantly broad quasielastic and inelastic peaks centered at hoj pa 0, 8 and 65 meV. Crystalline- 
electric-field (CEF) analysis satisfactorily reproduces the observed spectra, confirming their CEF 
origin. The broadness of the quasielastic peak indicates strong spin fiuctuations due to coupling 
between localized 4/ spins and conduction electrons in the paramagnetic phase. A prominent feature 
is suppression of the quasielastic fiuctuations, and concomitant growth of a sharp inelastic peak in a 
low energy region below Tn . This suggests dissociation of the conduction and localized 4/ electrons 
on ordering, and contrasts the presently observed incommensurate phase with spin-density-wave 
order frequently seen in heavy fermion compounds, such as Ce(Rui_a;Laa;)2Si2. 



I. INTRODUCTION 

The ternary rare-earth diborocarbidcs REB2C2 (RE: 
rare-earths) with heavy RE elements attract special at- 
tentions recently, as they exhibit intriguing successive or- 
dering at low temperatures. Exemplified by DyB2C2, 
the specific heat measurement detects two anomalies at 
Tc2 = 15.3 K and Tci = 24.7 K,i The higher temperature 
anomaly at Tci is attributed to ordering of quadrupole 
moments by various measurements such as the ultra- 
sonic vibration measurement^ and resonant X-ray difFrac- 
tioni^ii It has been suggested that the quadrupolar mo- 
ments originate from nearly degenerated two Kramers' 
doublets realized by pseudo-cubic local symmetry around 
the Dy'^^ ions. On the other hand, the lower tempera- 
ture anomaly at Tc2 is due to magnetic (dipole) order- 
ing. The magnetic structure established below Tc2 is 
non-trivial; a complicated multi-q ferrimagnetic struc- 
ture is realized with the four modulation vectors q = 
(1, 0, 0), (0, 1, 1/2), (0, 0, 0) and (0, 0, 1/2). Such a struc- 
ture cannot be stabilized by usual exchange-type inter- 
actions (including RKKY-type interactions), and strong 
influence by the quadrupolar order has been inferredii As 
exemplified above, interplay of dipolar and quadrupolar 
degrees of freedom is essential for the complex succes- 
sive ordering in the REB2C2 compounds with heavy RE 
elements. 

In contrast, REB2C2 with the light rare-earth element 
RE = Ce does not have quadrupolar degree of freedom; 
absence of the ground state degeneracy has been reported 
in the ultrasonic vibration measurement.— Nonetheless, 
CCB2C2 is of particular interest because influence of 
magnetic interactions to the complex successive ordering 
may be separately deduced from that of the quadrupo- 
lar interactions by comparing ordering behavior of non- 
quadrupolar CeB2C2 with those in the heavy RE sys- 



tems. Also, strong electron correlations originating from 
larger hybridization between 4/ and conduction electrons 
may bring about intriguing low-temperature properties. 

Earlier studies on CeB2C2 may be summarized as 
follows. Electronic specific heat coefficient 7 was esti- 
mated in the paramagnetic temperature range as 7 = 
98.7 mJ/mol K^i^ Thus, CeB2C2 may be classified as 
a heavy-fermion system with moderate mass enhance- 
ment, a typical consequence of strong electron correla- 
tions in the 4/ electron systems. A sharp anomaly at 
Tn = 7.3 K was observed in the specific heat and magne- 
tization measurements;^ indicating that a certain mag- 
netic order is established below this temperature. The 
specific heat anomaly is followed by a broad shoulder 
around T = 6.5 K, of which the origin is not elucidated. 
To date, only one attempt was made to solve the mag- 
netic structure; neutron diffraction experiment was per- 
formed using a single-crystal sample, and succeeded in 
observing a magnetic Bragg reflection at Q = {5, S, 6') 
with S w 0.161 and 6' w 0.100,1. On the basis of this 
observation, an incommcnsuratcly modulated structure 
was inferred in the ordered phase. However, because of 
insufficient quality of the crystal, it was the single peak 
that the diffraction experiment could detect, and con- 
sequently, the magnetic structure of the ordered phase 
could not be determined. Inelastic neutron scattering has 
also been carried out to determine the crystalline-electric- 
field (CEF) splitting of the 06^+ 4f levels*^ the inelastic 
spectrum was measured at 7 K for — 15 < fiij < 45 mcV. 
However, because of the limited energy range and tem- 
perature point, it hardly provides details of the spin fiuc- 
tuations and excitations in CeB2C2. 

Knowledge on the low-temperature magnetic structure 
is essential for the understanding of the magnetic order- 
ing in the CCB2C2, and its relation to REB2C2 with 
heavy RE elements. The CEF splitting and spin fluc- 
tuation spectrum are also crucial information to eluci- 
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date its magnetic properties. Hence, to address these 
issues we have undertaken low-temperature neutron pow- 
der diffraction and inelastic scattering experiments in 
the present work. In the powder diffraction we have 
succeeded in obtaining intensity data for a number of 
magnetic Bragg reflections, whereas inelastic spectra in 
a wide energy range up to 80 meV were collected at 
several temperatures below 100 K. These experiments 
enable us to unambiguously determine the magnetic 
structure in the ordered phase, as well as peak ener- 
gies in the inelastic spectrum; it will be shown that the 
lowest-temperature magnetic structure is a sinusoidally 
modulated structure with the modulation vector q = 
[0.167(3), 0.167(3), 0.114(3)], whereas the inelastic scat- 
tering detects broad quasielastic and two inelastic (CEF) 
peaks around huj « 0, 8 and 65 meV in the paramagnetic 
phase. The quasielastic fluctuations are suppressed be- 
low Tn, and a sharp inelastic peak emerges in the low 
energy region. 

II. EXPERIMENTAL DETAILS 

Polycrystalhne samples of Ce^^B2C2 (about 10 grams 
in total) were prepared using an arc furnace under an 



Ar gas atmosphere. Purity of the starting elements was 
99.9 % for Ce and C. To avoid strong neutron absorp- 
tion of natural boron, we used the isotope enriched ^^B 
(99.53 % enrichment). Structural quality of the resulting 
polycrystalhne samples was checked by the X-ray powder 
diffraction as well as the neutron diffraction. 



The polycrystalhne samples were crushed into powder 
and loaded in a thin Al sample cell for the neutron scat- 
tering experiments. The cell was then set to a '^He closed- 
cycle refrigerator. Neutron powder diffraction was per- 
formed using a powder-diffraction detector bank newly 
installed to the LAM-80ET inverted-geometry time-of- 
flight (TOF) spectrometer at the KENS spallation neu- 
tron source, KEK, Japan. The diffraction detector bank 
has 34 detectors covering the scattering angle range of 
138° < 29 < 162°. We used the incident neutrons in 
a wave-length range of 1.51 < A < 8.40 A. which cor- 
responds to the d-range of 0.8 < d < 4.3 A. Obtained 
powder diffractograms were analyzed using the home- 
made general-purpose Rietveld analysis code MSAS-TOFi^ 
To perform profile fitting of highly asymmetric peaks at 
large time-of-flights, a special profile function is used: 
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In the above expression, t' = fxoF — to, where to is a 
timc-of-flight at which a Bragg reflection appears. The 
two functions fcit') and fh{t') arc the Gaussian and 
Lorentzian functions defined by: 



foit') 
hit') 



exp(-t'V7^) 

7L 



- t'2) ■ 



(2) 



The profile for negative t' is a modified version of the 
pseudo-Voigt function; allowing independent values 
for 7G and 7l. On the other hand, to simulate a long 
tail due to the solid methane moderator at the KENS 
spallation source, a set of exponentially decaying func- 
tions is assumed for i' > as inferred from the Ikeda- 
Carpentcr equation]^ The parameters are assumed to 
be linearly to dependent in the present TOF range: 
V = Vi + V2to, 7g = 7G1 + 702*0, 7l = 7li + 7L2*o, 
ai = an -I- ai2to, 02 = 021 -I- a22to, "i = "ii + ^12*0, 
and a = (1 — rj) / (^/tt^q) + ri/{-K^i,). Details of the new 
diffraction bank at LAM-80ET and the analysis code will 
be published elsewhere i^ii^ 

Neutron inelastic scattering experiment was performed 



using the LAM-D inverted-geometry time-of-fiight spec- 
trometer, also installed at the KENS spallation neu- 
tron sourcCfi^ Final energy was fixed to 4.59 meV us- 
ing the pyrolytic graphite (PG) 002 reflections, whereas 
higher harmonic neutrons were eliminated by cooled Be 
filters. There are four analyzers at the scattering angles 
29 = ±35° and ±85°; we mainly show data taken with 
the lower angle detectors in this report, and higher an- 
gle data are used only for phonon subtraction purpose. 
Energy resolution was estimated as Ai? = 0.42 meV [full 
width at half maximum (FWHM)] at the elastic position 
using a vanadium standard. Background was subtracted 
from the raw data using a proper combination of empty- 
cell and absorber runs, and absorption effect was cor- 
rected using numerically calculated absorption factors. 



III. RESULTS AND DISCUSSION 

A. Crystal structure 

First of all, the validity of the newly installed detec- 
tor bank and the home-made Rietveld analysis code is 
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Parameters 


T = RT 


a (A) 


5.3948(3) 


c(A) 


3.8646(3) 


V (A3) 


112.47(1) 


XB 


0.3653(7) 


xc 


0.1614(7) 




5.4(3) 


Tpiso 


4.2(2) 


Diso 


3.5(2) 




0.0676 



TABLE I: Structural parameters for Ce^^B2C2 at the room 
temperature. The space group was assumed to be PA/mbm. 
The Ce atoms occupy the 2(a) sites at (0,0,0), whereas the B 
and C atoms occupy the 4(/i) sites at {xb,xb + 1/2, 1/2) and 
{xc,xc + 1/2, 1/2), respectively {Z = 2). 



confirmed by solving the room temperature structure of 
Ce^^B2C2. Figure 1 shows tlie powder diffraction pat- 
tern at the room temperature recorded without using 
the refrigerator. The pattern is analyzed assuming two 
phases; one is the tetragonal CeB2C2 structure with the 
space group PA/mhm determined by Onimaru et al.^ 
whereas the other is polycrystalline Al used as the sam- 
ple cell in the present experiment. For both the phases, 
the lattice parameters, atom positions (except for the Al 
phase), isotropic atom displacement parameters {B^^°), 
and preferred orientation parameters were optimized, in 
addition to the profile and background parameters. Re- 
sulting calculated diffraction profile is shown in the figure 
by the solid line. Difference between the observation and 
the calculation is also presented in the figure. Coinci- 
dence between the calculated and observed intensities is 
quite satisfactory. Obtained parameters are listed in Ta- 
ble I, which are in perfect agreement with the previous 
result This confirms the validity of the newly installed 
detector bank and the Rietveld analysis code. 

Figure 1(b) shows a powder diffractogram measured at 
T = 4.7 K < Tn. In this low temperature experiment, a 
serious background due to vacuum chamber walls of the 
refrigerator appears in a TOF region of 22.5 < iTOF < 
31 ms. Therefore, the data in the range were removed 
in the figure. It can be seen that the diffraction pattern 
below Tn is essentially the same as that at the room 
temperature. Hence, no change in the crystal structure 
is concluded in the present diffraction experiment. 



B. Magnetic structure 

Figure 2(a) shows the powder diffraction patterns at 
large TOF (30 < tTOF < 45 ms) in the low tempera- 
ture range (T < 8.0 K). Several additional reflections 
appear in the diffractogram below Tn. A huge Bragg 
peak at txOF = 33.9 ms is the superposition of the nu- 
clear 111 and 200 reflections. In the diffractogram, one 
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FIG. 1: (a) Powder diffraction patterns of Ce^^B2C2 at 
the room temperature. The calculated diffraction pattern is 
shown in the figure by the dashed line, whereas the difference 
between the calculated and observed intensities is depicted by 
the solid line, (b) Powder diffraction patterns of Ce^^B2C2 
at T = 4.7 K. The huge gap with no experimental data is 
due to the contamination of serious background due to the 
refrigerator chamber walls. 



can clearly see evolution of several new Bragg peaks be- 
low 7.0 K. They are two orders of magnitude weaker 
than the nuclear Bragg reflections. (Because of the weak 
intensity, they can hardly be recognized in the smaller 
TOF region.) Figure 2(b) shows temperature depen- 
dence of the integrated intensity of the Bragg peak ap- 
pearing at ixoF ~ 36.7 ms. On cooling, the intensity 
starts to increase at T sa 7 K, which is in good agree- 
ment with the macroscopically determined ordering tem- 
perature Tn = 7.3 K. This confirms the magnetic origin 
of the newly appearing Bragg peaks. We also find that 
these magnetic Bragg peaks can be indexed as satellite 
peaks of the nuclear reflections using the previously re- 
ported magnetic modulation vector, q = {5,5,5'), and 
its symmetrically equivalents. {5 and 5' are refined in 
the present study as described later.) Thus, the present 
powder diffraction result is consistent with the previous 
single-crystal study, and provides intensity information 
for more Bragg reflections that is mandatory for the spin 
structure analysis. 

Because of the limited TOF range and arbitrariness for 
magnetic reflection indexing due to powder averaging, it 
is not straightforward to find a spin structure model di- 
rectly from the observed diffraction pattern. Thus, to 
find possible spin-structure candidates, we use the mag- 
netic representation analysis introduced by Izyumov et 
al.^ In this method, structure candidates are given by 
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FIG. 2: (a) Neutron powder diffraction patterns at several 
temperatures between 8.0 K and 4.7 K. Magnetic Bragg re- 
flections clearly appear at the low temperatures. Dashed line 
stands for the result of the Rietveld fitting for the pattern 
at 4.7 K, whereas the solid line denotes the difference be- 
tween the observed and calculated intensities. Vertical lines 
shown in the lower part of the figure are the magnetic Bragg 
peak positions and intensities for magnetic structures given 
by the linear combinations of the irreducible representations 
lRn:m. Note that only IR2:1-|-IR2:2 can reproduce the ob- 
served diffraction pattern, (b) Temperature dependence of 
the integrated intensity for the (2 — 5,S,5') refiections, ap- 
pearing at tTOF ~ 36.7 ms. Solid line is a guide to the eyes. 



linear combinations of magnetic basis vectors of the irre- 
ducible representations in the paramagnetic phase. In the 
Landau theory of second order phase transition, a single 
irreducible representation may be selected as a symmetry 
of the ordered phase. In reality, it frequently happens 
that two or more irreducible representations are neces- 



sary to reproduce the symmetry of the ordered phase. 
Nonetheless, the number of the necessary representations 
is usually small. Thus, we may expect that the ordered 
phase in CeB2C2 may be given by a combination of a 
few irreducible representations. Here, wc try to find the 
magnetic structure model using the smallest number of 
the magnetic basis vectors. For this purpose, a repre- 
sentation analysis code, named MBASE, has been newly 
developed, which can calculate magnetic representation 
basis vectors for arbitrary fc-group and magnetic-ion po- 
sitions.— 

Assuming a single-q structure with multiple domains 
of the equivalent modulation vectors, direction of a spin 
(or total angular momentum) at the d-th site in the l-th 
unit cell may be generally written as: 



{Ji 



J 



l.d) 



[a<jcxp(-iq • Ri) + a;5exp(ig • i?/)], (3) 



where Ri denotes the position of the /-th unit cell. The 
polarization vector is given as a linear combination of 
the magnetic basis vectors: 



(4) 



where an^m.d denotes the basis vector of the irreducible 
representation IRn : m for the spin at the d-th site. The 
basis vectors for the q-domain and for the domains with 
symmetrically equivalent modulations (arms) arc listed 
in Table 11. In the present Rietveld analysis, the do- 
mains with the equivalent modulations are assumed to 
be equally populated. For the selection of the basis vec- 
tors, we note that there is very large anisotropy in the 
magnetic susceptibility; Xc is considerably smaller than 
Xa and Xiio-~ This indicates that spins most likely lie in 
the basal plane. Hence, we may use only the in-plane ba- 
sis vectors. In the bottom part of Fig. 2(a), reflection 
positions and intensities from the magnetic structures 
given by single or linear combinations of IRn : m are 
shown by the vertical thick solid lines. We find that a sin- 
gle irreducible representation cannot reproduce the peak 
positions; for instance IR1:1 definitely gives a peak at 
t-POF ~ 38.0 ms, which cannot be found in the observed 
diffraction pattern. Using two vectors with a constraint 
to the coefficient Cn,m ~ i'v/Z, wc find that the combina- 
tion IR2:1 -I- IR2:2 [i.e., a<j = (a2,i,d + a2,2,d)/V2], can 
satisfactorily reproduce the observed magnetic reflection 
positions. This combination corresponds to a sinusoidally 
modulated structure with spin polarization parallel to the 
[110] direction. The spin structure in the basal c-plane is 
schematically shown in Fig. 3. 

The profile fitting has been performed assuming the 
IR2:1-|-IR2:2 structure; the fitting parameters were the 
modulation vector q and amplitude J. Dashed line in 
the Fig. 2(a) shows the result of the profile fitting to the 
diffraction pattern taken at 4.7 K. The difference between 
the calculation and the observation is also shown in the 
figure by the solid line. Reasonable coincidence can be 
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IRn:m 




0-71,777,2 


q={S,5,S') 






1:1 


(1,0,0) 


(0,-e*,0) 


1:2 


(0,1,0) 


(-e*,0,0) 


1:3 


(0,0,1) 


(0,0,-e*) 


2:1 


(1,0,0) 


(0,e*,0) 


2:2 


(0, 1,0) 


(e*,0,0) 


2:3 


(0,0,1) 


(0,0, e*) 








1:1 


(1,0,0) 


(0,1,0) 


1:2 


(0,1,0) 


(1,0,0) 


1:3 


(0,0,1) 


(0,0,-1) 


2:1 


(1,0,0) 


(0,-1,0) 


2:2 


(0,1,0) 


(-1,0,0) 


2:3 


(0,0,1) 


(0,0,1) 


q = {S,~S,S') 






1:1 


(1,0,0) 


(0,1,0) 


1:2 


(0,1,0) 


(1,0,0) 


1:3 


(0,0,1) 


(0,0,-1) 


2:1 


(1,0,0) 


(0,-1,0) 


2:2 


(0,1,0) 


(-1,0,0) 


2:3 


(0,0,1) 


(0,0,1) 


q={-S,~S,5') 






1:1 


(1,0,0) 


(0,-6,0) 


1:2 


(0,1,0) 


(-6,0,0) 


1:3 


(0,0,1) 


(0,0,-6) 


2:1 


(1,0,0) 


(0,6,0) 


2:2 


(0,1,0) 


(6,0,0) 


2:3 


(0,0,1) 


(0,0,6) 



TABLE II: Magnetic representation basis vectors an^rn,d for 
the irreducible representations IRn : m. They are calculated 
for the four magnetic domains with the modulation vectors 
q = {5,S,5'),{—5,5,5'),{5, —5,5'), and (—5,-5,5'). e is de- 
fined as 6 = exp(27ri(5). The basis vectors for —q axe given by 
complex conjugates of the above vectors. The vector 
is for the (0,0,0) site, whereas a„,m,2 for the (i, |,0) site. 
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FIG. 3: Schematic drawing of the spin structure in the basal 
c-plane at T = 4.7 K determined in the present study. Only 
one domain with q = [5, 6, 5') is shown in the figure. 
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FIG. 4: Neutron inelastic scattering spectrum at T = 15.5 K. 
The open triangles and closed circles stand for the total scat- 
tering intensity and the phonon-subtracted magnetic inten- 
sity, respectively. Phonon contribution estimated from the 
T = 100 K data is shown by the solid line. 



C. Inelastic scattering 



1. Phonon subtraction procedure 



found in the figure between the calculated and observed 
intensities, despite the rather deficient statistics of the ex- 
perimental data. From the profile fitting. J is estimated 
as J = 1.8(2), corresponding to the maximum magnetic 
moment of gjJ^B = 1-5(2) Mb at 4.7 K. The q- vector 
is refined as q = [0.167(3), 0.167(3), 0.114(3)]. It may be 
noteworthy that this q-vcctor is quite close to the com- 
mensurate value of (1/6, 1/6, 1/9), although this cannot 
be concluded in the present powder diffraction measure- 
ments; a single-crystal experiment is highly desired to 
conclude this issue. 



Next, to investigate the spin fluctuations and excita- 
tions, we have performed the inelastic scattering experi- 
ment using the same powder sample. To reliably obtain 
magnetic scattering intensity from the raw inelastic spec- 
trum, a phonon contribution has to be carefully removed. 
Hence, we first make an estimation of the phonon con- 
tribution using high temperature data at T = 100 K. 
For the rare-earth compounds in a paramagnetic tem- 
perature range, magnetic scattering may originate from 
local (single-site) transitions, and thus, Q-dependence of 
the magnetic scattering may be dominantly given by the 
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FIG. 5: Magnetic scattering spectra at T = 52.05, 32.3, 23.25, 
and 15.5 K (from top to bottom). The solid lines stand for 
the fitting to the model scattering function Eq. (|8} , combined 
with a delta-function like elastic peak. See text for details. 
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FIG. 6: Magnetic scattering spectra at T = 52.05, 32.3, 23.25, 
and 15.5 K (from top to bottom). The solid lines stand for 
the fitting to the pseudo-Voigt quasielastic signal, and a delta- 
function hke elastic peak. See text for details. 



magnetic form factor. On the other hand, the phonon 
scattering may approximately exhibit dependence. 
Thus, we assume the following Q dependence for the total 
(ie., magnetic -I- phonon) scattering fimction: 

5tot(Q, fi^) = Anag(fi^)[./mag(Q)]' + V(fiw)g^ (5) 

where A-^g_^(hiS) and Aph(?iw) stand for the Q inde- 
pendent parts of the magnetic and phonon scattering, 
whereas /mag(Q) for the magnetic form factor of the 
Ce^+ ionsiii By comparing the inelastic spectra at T = 
100 K measured with the two different scattering angles 
(26' = 35° and 85°), we estimate the phonon contribution 
A-p\Jyf}w)Q'^ . Then, the phonon contribution at low tem- 
peratures is obtained using the [1 — exp(— ftw/fceT)]^^ 
dependence. Figure 4 exemplifies this phonon subtrac- 
tion procedure for the representative data at T = 15.5 K. 
The solid line stands for the estimated phonon intensity, 
whereas phonon-subtracted magnetic intensity is shown 
by the filled circles. In the low energy region, the phonon 
intensity is considerably small, whereas the phonon in- 
tensity becomes dominant in the higher energy region 



(hw > 30 meV). This large phonon contamination re- 
duces statistical accuracy of the estimated magnetic in- 
tensity, however, one finds that there definitely exists fi- 
nite magnetic intensity for huj > 50 meV. 



2. Inelastic spectrum in the paramagnetic phase (T > Tisi) 

The phonon-subtracted magnetic scattering spectra at 
four representative temperatures T = 52.05,32.3,23.35, 
and 15.5 K in the paramagnetic phase are shown in 
Figs. 5 and 6. Three peaks are recognized in the spec- 
tra at all the temperatures, in addition to the sharp 5- 
function like elastic peak: (i) quasielastic component cen- 
tered aXhw = Q (Fig. 6); (ii) inelastic peak at ftw sa 8 mcV 
(Fig. 5); (iii) weak inelastic peak at fiij w 65 meV 
(Fig. 5) . Since those spectra were measured in the param- 
agnetic phase, the peaks are most likely due to transitions 
between the CEF splitting levels. Hence, the observed 
spectra are analyzed using the following CEF Hamilto- 
nian derived for 00"^+ in the C4/1 site symmetry of the 
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2(a) site in CeB2C2: 

^^CEF = B20O20 + B40O40 + -B44O44 + i?4_404_4, (6) 

where 0„,„ stands for the Stevens operators Under 
above CEF, J = 5/2 multiplet of Ce^+ spHts into three 
Kramers' doublets. The transition strengths between the 
CEF spUtting levels are given as follows: 

2e-^-A-BT|(^|jo|^)|2 
Km = ^ ^ (m t n) 

= ^ \^rj, , (otherwise) (7) 

I 



where |n) and \m) are wave functions for the initial and 
final states of the CEF splitting levels (see Fig. 7 for 
numbering of the states) , and Z is the partition function. 
The scattering function from a powder sample may be 
given by a sum of spectral weights of the CEF transitions: 



SiQ,hw) 



incl 



: 3,7 /mag (Q) 



Nhuj 



1 — cxp(— /la^/fceT) 



In the above equation, we assume a Lorentzian-type peak profile for the inelastic CEF excitations: 



1 



1 



I 



(8) 



(9) 



For the quasiclastic peak shape, we assume a pseudo- 
Voigt function, which is a reasonable approximation of 
Lorentzian function convoluted by a Gaussian-shaped in- 
strumental resolution functionJ^iii For the pseudo-Voigt 
function, the width of the unconvolutcd Lorentzian is 
denoted by r„„ (FWHM). We note that the lowest or- 
der coefficient ^20 can be determined precisely using a 
single crystal magnetization measurement; -B20 is esti- 
mated as 6.34 K in the earlier worki^ Thus, we fix ^20, 
and try to find the optimum values for i?4o, -B44 and 
-84-4 that reproduce all the inelastic spectra in the para- 
magnetic phase simultaneously. In the fitting, we as- 
sume that temperature dependence of the Hamiltonian 
parameters -820 j ^40 j ^44 j and -64-4 is negligible in the 
present temperature range. On the other hand, most of 
the width parameters r„m are assumed as temperature 
dependent; only Fia and r23 have to be fixed to the em- 
pirical value 12 meV because of the insufficient statistics 
in the high energy regions. The resulting optimum pa- 
rameters for the CEF Hamiltonian arc -640 = —1.84(3) K, 
B44 = 10.1(2) K and ^4-4 = -2.8(6) K. Calculated spec- 
tra using these parameters are shown in Fig. 5 by the 
solid lines. Reasonable coincidence can be seen between 
the observed and calculated intensities at all the tem- 
peratures. The fact that the temperature dependence is 
well reproduced confirms the CEF (magnetic) origin of 
the three peaks. It should be noted here that parame- 
ter -64-4 could not be estimated reliably; its uncertainty 
range is considerably larger than those of other param- 
eters. This parameter exists in the CEF Hamiltonian 
because of the broken four fold symmetry due to the B- 



C ordering. Since B and C have relatively similar elec- 
tronegativity, the symmetry breaking may possibly be 
moderate, and thus we might infer that i?4_4 may be 
irrelevant. 

Earlier neutron inelastic scattering experiment pro- 
vides considerably different Hamiltonian parameters; 
S40 = 0.024(2) mcV (= 0.28 K) and B44 = 0.356(1) meV 
(=: 4.13 K)4 The discrepancy is due to the assignment 
of the highest inelastic peak; the previous study assigned 
the highest energy peak to a very broad hump found at 
hw sa 23 meV in their spectrum. The hump is completely 
absent in the presently observed inelastic spectrum. It 
may be noteworthy that a weak peak in phonon density 
of states of the elemental aluminum, which is commonly 
used for sample cells in inelastic neutron scattering, ex- 
ists in this energy range, and thus this may become a 
source of uncertainty. On the other hand, we clearly see 
the highest energy peak at ftw ~ 65 meV, which is out of 
the observation energy range of the previous study. We 
also note that earlier macroscopic study provides rough 
estimate of the CEF levels at 102 K and 1110 sup- 
porting our result. 

Among the three peaks, the quasiclastic peak is of par- 
ticular interest in the Ce compounds, since it reflects the 
low-energy scattering process between the ground state 
doublets and the conduction electrons. Hence, we pa- 
rameterize the quasielastic-peak width as a function of 
temperature. To obtain the width precisely, the fitting is 
performed in the limited energy range —3 < fio; < 6 meV 
using the pseudo-Voigt function with the unconvolutcd 
Lorentzian width Fqoi as fitting parameter. The fitting 
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FIG. 7: Schematic illustration of the CEF splitting scheme in 
the CeB2C2. The energy separations of the Kramers' doublets 
are determined in the present study. See text for details. 
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FIG. 8: Temperature dependence of the quasielastic peak 
width Fqpi (FWHM) in the paramagnetic temperature range. 
The solid line is a guide to the eyes. 



result is depicted by the solid lines in Fig. 6. Obtained 
temperature variation of the quasielastic peak width Fqci 
is shown in Fig. 8. The quasielastic peak becomes nar- 
rower on cooling, as generally seen in heavy-fermion sys- 
tems. It may be noteworthy that the quasielastic width 
is noticeably large even at T w 12 K, just above Tn. 
This indicates that quasielastic fluctuations due to the 
scattering by the conduction electrons are still dominant 
in vicinity of the ordering temperature. 
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FIG. 9: (a) Neutron inelastic spectra at T = 12.27 K > Tn 
and T = 4.77 K < Tn. Solid line is the calculated spectrum 
assuming the CEF Hamiltonian Eq. ((T]), whereas the dashed 
line for T = 4.77 K is the result of the scattering-function 
calculation assuming an additional Zeeman term in the CEF 
Hamiltonian. (b) Excitation spectra in the low energy region. 
Solid lines in this plot are guides to the eyes. 



3. Inelastic spectrum in the ordered phase (T < Tn ) 

As seen in the previous section, the temperature de- 
pendence of inelastic spectrum is found to be moderate 
in the paramagnetic phase. However, it shows drastic 
change across Tn. Representative inelastic spectrum be- 
low Tn is shown in Fig. 9, in comparison with the param- 
agnetic spectrum at T = 12.27 K. The broad quasielastic 
signal around fuo — Q disappears below Tn, and a new 
inelastic peak develops a,t fku = 2.1 meV, instead. Note 
that the newly appearing 2.1 meV peak is quite sharper 
than that of the quasielastic peak at T > Tn- Thus, 
the low-energy spin fluctuations in the ground state dou- 
blet are strongly suppressed on ordering. This strongly 
suggests dissociation of 4/ electrons from the conduction 
electrons, i.e. formation of the localized moment in the 
ordered phase. 

On the origin of the sharp inelastic peak, several pos- 
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sibility may be anticipated. The simplest possibility may 
be splitting of the ground state doublet by the inter- 
nal (exchange) molecular field appearing in the ordered 
phase. To check this possibility, we calculated the single- 
site CEF excitation spectrum under the internal field by 
introducing the Zeeman term iJzocman = QjI^bJ ■ -ffint 
in the CEF Hamiltonian Eq. ([7]). Direction of the inter- 
nal field Hint is parallel to the spin direction [110]. The 
calculated inelastic spectrum assuming Hiat 17 T is 
shown in Fig. 9(a) by the dotted line. The peak position 
is reproduced by introducing the internal field, however, 
relative spectral weights of the peaks are apparently in- 
consistent with the observed spectrum. Since in the si- 
nusoidally modulated structure the molecular field is not 
uniform, another possibility may be the doublet split- 
ting due to the distributed (non-uniform) internal fields. 
However, the distribution will apparently introduce peak 
broadening, and thus will not reproduce the sharpness 
of the 2.1 meV peak. As above, the single-site CEF ori- 
gin is unlikely for the sharp inelastic peak. As another 
possibility, we may thus speculate that the inelastic peak 
plausibly stems from a certain collective excitation of the 
interacting Ising-like spins formed by the ground state 
doublets. To pursue the origin of the 2.1 meV inelastic 
peak, further study on its Q-dependence using a single 
crystal is necessary. Such a single-crystal inelastic scat- 
tering experiment is in progress. 

It should be reminded that the long-period sinusoidally 
modulated structure has a number of spins that have 
strongly reduced average spin magnitudes (S) « S (see 
Fig. 3). In a localized spin system, such reduced spins are 
realized by thermal fluctuations, and thus becomes un- 
stable for further lowering temperature (22. On the other 
hand, for an incommensurate phase in hcavy-fcrmion sys- 
tems, such reduced spins are formed by quantum fluctu- 
ations; coupling with conduction electrons enables a for- 
mation of spin-density-wave-type (SDW) order, as typi- 
cally seen in Ce2:Lai_:rRu2Si2^ or Ce(Rui_i;Rh2;)2Si2'^ 
In both the cases, either thermal or quantum spin fluc- 
tuations should remain in the sinusoidally modulated 
phase. In contrast to the above understanding of the si- 
nusoidally modulated phase, quasielastic fluctuations are 
completely suppressed in CeB2C2. This is a very unique 
characteristic of the sinusoidally modulated phase in the 
CeB2C2, and further experimental as well as theoretical 
study is highly desired to clarify this issue. 

Finally, we compare the ordered spin structure of 
CeB2C2 with those in the heavy RE systems. Incom- 
mensurately modulate structures have been frequently 
observed in the heavy RE systems, such as IIoB2C2r^i^ 
TbB2C2^'^ and ErB2C2-21 In H0B2C2, incommensu- 
rately modulated spin structure is realized blow Td = 
5.9 K with the two modulation vectors q = (1,0,0) and 
(1 ± (5, S, S') where S = 0.112 and 6' = 0.04. Quadrupolar 
order is established at lower temperatures T < Tc2 ~ 
5.0 K, where spin structure becomes commensurate with 
the four modulation vectors q = (1,0,0), (0,1,1/2), 
(0,0,0) and (0,0,1/2). Note that the spin modulation 



vectors in the quadrupolar ordered phase are identical to 
those in DyB2C2. TbB2C2 also shows similar multi-q in- 
commensurate structure. This compound does not show 
quadrupolar ordering in the zero external field, however, 
quadrupolar ordering is known to take place in a very 
low external magnetic field of 1 T, indicating that a likely 
situation is realized for the formation of the quadrupolar 
moment in the CEF ground state. In contrast, ErB2C2 
exhibits the singlc-q sinusoidally modulated structure be- 
low Tn = 15.9 K with the incommensurate modulation 
vector q = {1 + 5,6,0) {S = 0.112). It is followed by 
a lower-temperature lock-in transition to q = (1,0,0) 
at Tt = 13.0 K. Such a lock-in transition is typical for 
the linearly polarized incommensurate phase in localized 
spin systems. No quadrupolar ordering behavior was re- 
ported in this system. In the present study, we have also 
observed a single-q sinusoidal phase in the CeB2C2 with 
no quadrupolar degree of freedom. Therefore, this re- 
sult, in addition to the above observations for the other 
RE systems, supports the claim that the magnetic sec- 
tor of the REB2C2 compounds has a tendency to form a 
long-period single-q structure, and it is the effect of the 
quadrupolar degree of freedom that realizes the complex 
multi-<7 structures. 



IV. CONCLUSIONS 

We have performed neutron powder diffraction and in- 
elastic scattering experiments on the cerium diborocar- 
bide Ce^^B2C2. In the powder diffraction study, we have 
clearly detected the magnetic reflections below the or- 
dering temperature Tn = 7.3 K. With the aide of the 
magnetic representation analysis, we have determined 
the magnetic structure in the ordered phase as the sinu- 
soidally modulated structure with the modulation vector 
q = [0.167(3), 0.167(3), 0.114(3)]. The inelastic study de- 
tects three magnetic peaks in the paramagnetic phase; 
the quasielastic signal around Huj = and broad in- 
elastic peaks at hu! ^ 8 and 65 meV. Upon cooling to 
T < Tn, drastic change has been seen in the low energy 
spin fluctuation spectrum; a sharp inelastic peak develops 
at hu! = 2.1 meV. This indicates that the spin fluctua- 
tions due to hybridization between the 4/ and conduction 
electrons in the paramagnetic temperatures are strongly 
suppressed in the ordered phase. 
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